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We present a first A'^f = 2 lattice estimate of the hadronic coupling gi2 which parametrises the 
strong decay of a radially excited B* meson into the ground state B meson at zero recoil. We work 
in the static limit of Heavy Quark Effective Theory (HQET) and solve a Generalised Eigenvalue 
Problem (GEVP), which is necessary for the extraction of excited state properties. After an extrap- 
olation to the continuum limit and a check of the pion mass dependence, we obtain gi2 = —0.17(4). 

PACS numbers: 12.38.Gc, 13.20.He. 



55^. I. INTRODUCTION 

Questions have been raised recently on the poor handling of excited states in the analyses of experimental 

data and their comparison with theoretical predictions. For instance, it has been advocated that the ~ 3a 

^ . discrepancy observed between V^'^ff'^' and F^™'^' may be reduced if the transition B ^ D' were large. This 



cb ch 

^^ \ attractive hypothesis implies a suppression of the B — ^ D^*> hadronic form factors, as a study in the OPE 

formalism suggests jlj] . On the other hand, it has been argued that the light-cone sum rule determination 
of the go'D-n coupling, which parametrises the D* — >• D-k decay, likely fails to reproduce the experimental 
measurement unless one explicitly includes the contribution from the first radial excited D^*^ state on the 
^J ' hadronic side of the three-point Borel sum rule [2| . Comparison with sum rules is of particular importance 

(>^ . because the heavy mass dependence of gq = 2^hn"m''f deduced from recent lattice simulations [j-Q and 

experiment [9| seems much weaker than expected from analytical methods 10], as shown in Figure [Ij 

Techniques have been developed to study excited states of mesons using lattice QCD [11], especially to 
extract the spectrum |12h15| |. Similar techniques can now be applied to three-point correlation functions 
to perhaps illuminate the phenomenological issues discussed above. In this letter we will report on the 
lattice computation of 512 = {B* \Ai\B) in the static limit of HQET, where Ai is the axial vector bilinear 
of light quarks and B* is polarised along the ith direction. As a by-product of our work, we will also 
report on the computation of gn = {B*\Ai\B) and 322 = {B* \Ai\B'). 

The Heavy Quark Symmetry of leading order HQET is well suited for our qualitative study. As the 
spectra of excited B and B* mesons are degenerate, it is enough to solve a single Generalized Eigenvalue 
Problem (GEVP) while degrees of freedom ~ m\j, that are somehow irrelevant for the dynamics of the 
cloud of light quarks and gluons that governs the process we examine, are integrated out. The plan of the 
letter is the following: in Sec. |ll]we describe our approach while in Sec. IIIII we present our lattice set-up 
and discuss results before concluding in Sec. IIVI 
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FIG. 1: Experimental measurement [y|, lattice computations [3431 and sum rules estimates [Ifll oi Qc, Qb and g = g^a. 
We have performed a weighted average of recent g lattice results at Nf = 2 with respect to the error quoted in 



II. EXTRACTION OF {B* \Ai\B) 



The transition amplitude of interest is parametrised by 

c(A) . 



q^ V r 

[pb+Pb"Y 



g(A) .0 
+ A,{q^)- '^ 



m-B + nig,' 



(1) 



with q = p' —p. In the zero recoil kinematic configuration where p ~ p =0, one has (?niax — {rrig,' — tub)'' 



so that 



{B*'{p',e,)\A^\B{p)) = {mB + mg.,)A,{ql,^)e^^^\ 



(2) 



At that stage it is useful to introduce the HQET normalisation of states: \H) — y^2mH\H)iiQET, with 
{H{p)\H{p'))^2E{p)6^{p-p): 



2y/mBmg,' 



(3) 



In the static limit we are left with {B* (p', eA)|-4'|B(p))HQET = ^i('?max)^'''^'''- Choosing the quantization 
axis along the z direction and the polarisation vector e^(0) = n i with the metric (+,-,-,-), we get 

finally Ai{q'^-^^^) = {B* (p', eo)|y^'^|i?(p))HQET- Of course, extracting gu = g and 522 is similar, except that 
the relevant axial form factors are defined at q^ = 0. 

GEVP methods [16l - [l8| are a very efficient tool to study excited states on the lattice. We consider 
N X N matrices of two-point correlation functions together with the corresponding matrices of three- 
point correlation functions C^^-^^^^) = {Op{t)Op\o)) and C^fit,t,) = {O[(t,)Or{t)O]{0)), where i,j 
represent different wave functions and Dirac structures with quantum numbers generically denoted (h). 
More explicitly, the Oi are interpolating fields of pseudoscalar static-light mesons, the O'^ interpolating 
fields of vector static-light mesons and Or the axial vector light-light bilinear of quarks. 



In HQET the spectral decomposition reads cfp^^^(i) = E„ Cf '''Vi'/'^e"-^"*, vif'' = 

{Mn \Oi |0). The purpose of solving GEVP is to construct quantities which tend toward the desired 
excited state properties asymptotically in time. In practice we solve 

3 3 

We will use two ratio methods, GEVP and sGEVP, to extract the matrix element Mm„ = (M„ ^ |OrJMm '). 
Those ratios converge quickly as the contribution of higher excited states is strongly suppressed jl9|^ and 
read: 

™" ' v/('^("^ {is - t, to),Ci^^ {ts - t)vM {u - i, to)) (^'(") {t, io), C'(^) (t)«(™) (t, io)) 

i 

^sGEVP.^^ ^ 5 r (i;("')ft,to),[i^'""(t,to)/A(")(t,to)-J^— fa,to)]z;(")(t,to)) 

" * [V(«(™'(i,io),i?™"(t,io)«(™'(i,io))V(«("Hi,io),C(2)(io)«(«)(i,io)) 

~ Mmn+0(Ate °), (6) 

*i 

In the appendix, we have calculated the time dependence of the corrections in R^^^^ (i) to first order in 
e, where 

AT 00 

n=l n=N+l 

nsGEVP _ 7i,r , ,nsGEVP,l 

We have found that for n > m the dominant contribution to eR^^^^'^ is ie~(^"+i~^"^* and for n < m 
the leading contribution is in e-iEN+i-E^)t^ 

The global phase is fixed by imposing the positivity of the 'decay constant' /^,^(h) = (Af„ ^|O{^|0) — 



^{v(^Ht,to).C(^Ht)v\"\t,to)) 



-, where L refers to some local interpolating field. 



III. LATTICE RESULTS 



We have performed measurements on a subset of the Nf = 2 CLS lattice ensembles, which employ the 
plaquette gauge action and non-perturbatively 0{a) improved Wilson-Clover fermions. The parameters 
of the ensembles used in this work are collected in Table ID Three lattice spacings (0.05 fm < a < 0.08 fm) 



We give in the Appendix a hint of the proof of the t behaviour of R^^^^ (t) , as it was not discussed in detail in [l9ll . 



0.28 1 



0.24 - 



f^ 



^ 



w" 0.2 ■ 



0.16 ■ 



0.12 - 



12 



t/a 



FIG. 2: Plateau of E2 - Ei for the CLS ensemble E5. 
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FIG. 3: Dependence of bare gi2 on the size of the GEVP (left) and on the radius of wave functions (right) for the 
CLS ensemble E5. 



are considered with pion niasses in the range [310 MeV, 440 MeV] . The static-hght correlation functions 
employ the 'HYP2' discretization of the static quark action [20, l21| and stochastically estimated all-to-all 
light quark propagators with full time dilution [24I . A single fully time-diluted stochastic source has been 
used on each gauge configuration, except for the ensemble E5 where we have four stochastic sources for 
each gauge configuration. We use interpolating fields for static-light mesons of the form [23| 



where kg = 0.1, r^ = 2a^ KQRi < 0.6fm, and A is a gauge covariant Laplacian made of 3 times APE- 
blocked links ^. 

In order to reduce the statistical uncertainty in ratio (jH]), we have taken the asymptotic value of the 
energy splittings S™" = En — Em- We have shown in Figure [5] an example plateau for S^. In addition 
we have set ts to 2t in ([5]). We have solved both 3x3 and 4x4 GEVP systems and checked the stability 
of the results when the local operator is included, as shown in Figure [H To check the dependence on ig, 
to which the contribution from higher excited states is sensitive, we have both fixed it at a small value 
(typically, 2a) and let it vary as i — a. 

Though the uncertainty is a bit larger, we have confirmed the finding by [l^l that using sGEVP ((H) 
seems beneficial compared to the standard GEVP approach ([5]) to more strongly suppress contamination 
from higher excited states in the hadronic matrix element we measure. As illustrated in Figure 21 plateaux 
obtained from the GEVP and sGEVP are compatible: -0.25(1) for GEVP and -0.23(2) for sGEVP, with 
one additional point in the plateau of the sGEVP. Therefore, in the following we give results using the 
sGEVP only. 

After applying a non-perturbative procedure to renormalise the axial light-light current |25l . l26j . we 
are ready to extrapolate to the continuum limit. Inspired by Heavy Meson Chiral Perturbation Theory 
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FIG. 4: Plateaus of bare gi2 extracted by GEVP (left) and sGEVP (right) for the CLS ensemble E5. 

at leading order [27, 28] and due to the 0{a) improvement of the three-point correlation functions (the 
improved part of the axial current, acAdiV, is absent at zero momentum), we apply two fit forms: 



512 = Co + (0/0/3=5.3)^^1 + (m^/m^)^C2 , 
912 = C^ + (a/a/3=5.3)^Ci • 



(7) 
(8) 



We show in Figure [5] the continuum extrapolation ([7]) of gi2. We observe quite large cut-off effects (~ 30% 
at ;S = 5.3), it is thus crucial to have several lattice spacings. We obtain finally, using (O as the best 
estimate of the central value, 



gi2 = -0.17(3)(2), 



(9) 



where the first error is statistical, and the second error corresponds to the chiral uncertainty that we 
evaluate from the discrepancy between ^ and ^. We collect in Table HIl the value of 312 at each lattice 
point and at the physical point as well as the fit parameters for ([7]) and ([8]). 

In simulations with light dynamical quarks, the onset of multi-hadron thresholds due to the emission 
of pions must be considered when examining excited B meson properties. Such thresholds significantly 
complicate the extraction of hadron-to-hadron matrix elements from the two- and three-point correlation 
functions considered here. However with the L < 3 fm volumes in this work, the P-wave decay B* (0) — > 
B(p)'k(—p) is kinematically forbidden. The S-wave decay B* — >■ B^t: is potentially more dangerous. 
Examining the mass splittings E12 in Table III, we notice that 630 MeV < S12 ^ 710 MeV. If we assume 
that 400 MeV < ms. - tub < 500 MeV in the pion mass range [310 MeV, 440 MeV], (as has been found in 
a recent lattice study of the static light meson spectrum |29|). we conclude that our analysis is safe from 
these threshold effects. Moreover the bare couplings gi2 we obtain are similar to the quenched result of 
Ref. [H. 

We show in Figure [6] a typical plateau of the bare coupling gn and the extrapolation to the continuum 
and chiral limit. That extrapolation is smooth, with a negligible dependence on TTi^r, and we obtain from 
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TABLE I: Parameters of the simulations. 





312 


A5 
E5 
F6 
N6 


-0.245(29) 
-0.186(8) 
-0.207(15) 
-0.181(12) 


physical point 


-0.173(28)(18) 





fit (III) 


fit (18|) 


Co 
Ci 
Ca 


-0.178(29) 
-0.063(32) 
0.0053(29) 


-0.155(26) 
-0.040(29) 



TABLE II: Value of gi2 at the lattice points and at the physical point (left) as well as the fit parameters of Eq. ((TJ 
and (O (right). 





aEi2 


A5 
E5 
F6 
N6 


0.255(8) 
0.222(8) 
0.216(12) 
0.173(7) 



TABLE III: Mass splitting S12 in lattice units. The error we quote is the discrepancy between plateaux that we 
extract for different time ranges {[fmin, imax], [tmin ± 0.2ro, imax ± 0.2ro]}. 

the fit form ([7]) gu = 0.52(2), in excellent agreement with a computation by the ALPHA Collaboration 
focused on that quantity Q. We have added an error of 2% due to higher excited states which is estimated 
from plateaux at early times with a range ending at ^ rg. Following the same strategy, we show in 
Figure [7] a typical plateau of the bare coupling 322 and the extrapolation to the continuum and chiral 
limit, once again quite smooth, with an almost absent dependence on the sea quark mass. We obtain 
from the fit form ([5]) (722 = 0.38(4). Remarkably, the "diagonal" couplings gu and (722 are significantly 
larger than the off-diagonal one gi2 ■ This suggests that neglecting the contribution from B' mesons to the 
three-point light-cone sum rule used to obtain gB'Bir introduces uncontrolled systematics. Note that the 
decay constant fg,i itself is large compared to Jb [30,L31]. For completeness we have collected in Table HVl 
the value of 51 1 and (722 at each lattice point and at the physical point and the fit parameters of ([7]) and 
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FIG. 5: Continuum and chiral extrapolation of g\2 
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A5 
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F6 
N6 


0.541(5) 
0.535(8) 
0.528(4) 
0.532(6) 
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0.474(26) 
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g22- fit ISJI 



0.385(24) 
0.076(26) 



TABLE IV: Value of gn and 322 at tlie lattice points and at tfie pliysical point (left) and fit parameters of eq. ((TJ 
and ((HI) (rigfit). The third error on gn is an estimate of the effects of higher excited states. 
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FIG. 6: Plateau of bare gw for the CLS ensemble E5 (left) and its extrapolation to the continuum and chiral limit 
(right). 



IV. CONCLUSION 



We have performed a first estimate of the axial form factor A\{q^^^ = gi2 parametrising at zero recoil 
the decay B* — > _B in the static limit of HQET from Nf = 2 lattice simulations. Assuming the positivity of 
decay constants fs and fg,' , we have obtained a negative value for this form factor. It is almost three times 
smaller than the gn coupling: 512 = —0.17(4) while gu — 0.52(2). Moreover we find 1722 = 0.38(4), which is 
not strongly suppressed with respect to gu. Our work is a first hint of confirmation of the statement made 
in Ref. [2| to explain the small value of go'D-K computed analytically when compared to experiment. This 
computation using light-cone Borel sum rules may have been too naive. Following Ref. [32|, a next step in 
our general study of excited static-light meson states would be the measurement of Ai (0) by computing 
the distribution in r of the axial density fAir) = {B*' \i[)i"f^"f^il:i{r)\B) and Ai(0) = 47r j^ r'^ fA{r)e^^''^ dr. 
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Appendix 

In this section we discuss the time dependence of -R^^^^ (HI) • To simplify notation, we have fixed a to 
1. We have followed the strategy of Ref. JJi] to treat in perturbation theory the full GEVP, with an exact 
computation of the N lowest states: 

N 00 

Clf (i) = (O.(t)O,(0)) = C^^^it) + ^-^^{i) ^ Y. ^-''''Vn.V'n, + E e-^'^V^V'n, , 

n=l ra=JV+l 

C(2)(i)i,„(i,to)=A„(t,io)C(2)(ioK(t,to), 
w„(t,io)-wi°'(^,io) + CTi'Hi,io), 
A„(i,io)-Al")(i,to)+eAW(t,io). 



Vectors are normalised such that 






where /0„ = e ^"*" . Introducing the dual vectors u„ defined by X]n=i '^ni'4'mi = ^mn Vn < iV, we note 
that 

At first order in e, we have 



V. 



/^ ^m A'm (0) _ , (0) Z^ 



L^nm. *-^^ 



where A„ = C(2.i)(i) - X^ {t,to)C^^^^Hto)- With c„,™j = (M„,Vi )("«>;) we get: 



An {t,to) 



An (ij^o) /' 



/ ^ ^n,n^l^ 



-{Ei-E„)to 



I _ g-(Si-£;,.)(t-to) 



>Af 



e anm(i,^o) — ~ ^ Cn^mJ- 



1 „ g-(£;,--E„)(t-to) 



-(£;,-B„.)*o 



;>A' 



= -(£;,„-B„)(t-io) 



Finally the normalisation conditions read 
We are ready to develop ^ to first order in e: 



M-f;: = a, 



{Vmjt, to), [K{t, to)/Xn{t, to) - K{to, tp)] Vnjt, tp)) ^S{to,to) 



[{Vrn{t, to), C(2) {to)Vmit, to)){Vnit, to), C(2) (io)^n(i, to))] 



1/2 



e 2 






A<$;f„^'° - ft 



y^m-! 



K{t,to)/X^°\t,to)-K{to,to) 



[(u„,C(2)(io)u„)(^*„,C(2)(to)^t„)] 



1/2 



,^E(to,to) 



With 



{um,K{t,to)un) = ^e-^(*-*i)(w™,C(3)(t,ii)u„) = /i^n te-^"*, (u„,C(2)(to)u„> 



Pn ^ e 



^E„to 



we have at leading order 
The subleading order reads 



A^eff,s,0 _ L 



eM^±^' = edtJ^Ta- 



,(0)/ 



Ti = _ A (Mo) (.^(Mo),J^(t,M^n (Mo)) ^M.(.„,„) ^ -^#^(x'L"'(Mo),if(Mo).i"'(Mo)>. 

(Ak^(i,io))2 (PnP™)l/2 Xr{t,to) 



E„t 



The first subleading contribution is given by 



-t 1 — "mn ''Co) <-«,«, Af+1 "mn ^ (-C 

An {t,to) 



1-e 



^A7v + i,Ti(i — io) 



Defining the discrete derivative dtA~ A{t + 1)~ A{t), and taking at the end of the computation tp — t — 1, 
we get 



dtTi ~ C„,n,7V+l hmn (l " e-^" + i'") X (i + 1 + 

The second subleading contribution reads 

„(i) 



gAjv + i,,i\ g-Ajv + i,„t 



T, 



(i, to), [K{t, to)/A.L"^ (t, to) - K{to, tp)] v^^^ {t, to)) ^ 



(PnPm)^/^ 



-S(to,to) 



^ a™p («^*"(t,to), [i^(t,to)e^"* - A'(to,to)e^"*''] «i°^(t,to)) 



10 



With some algebra, we deduce 

ti rs 

^ y^g-(t-ti)(iJ„--E™)/j ^g--Ep(t-ii)g--E„tig£;„t 



J2hp 



-(£p-B„)ti 



and 



ti=to+l 



-E^)ti 



Finally, 






o{t,to) ^ hp 



-(Ep-E^)ti 



ti=to + l 






ti=to+l 



Setting to = t — 1, the first term reads 

- 1 _ ^-2(Em + i-E^) ^ _ ^-(En + 1-E^) 



E 

p^m 



'-■va.p.N+l^ ^ 



X Ve-^^'-^"^* 



_g-(B„+i-B™)t y- 

p^m 



and the second term reads 






I _ g-2(i3p-_E„) ^ _ g-(Bp-B™) 

^ g— 2{£'iv + i — -Etti) ]^ ^ — {ENJ^i—Em) 



I _ g-2(Bp-B„) I _ g-{Ep-E„ 



^a„,(i + l,to) Ve-(^''-^'")(*+^^ 



p^jn 



■E^ 

p^m 



1 „— 2(£'jv + i — -Etti) 

-(i3« + i--Ep)(t-l).!__f „ ^^,^^ -(£p-i3„)(t+l) 

I _ g-2(Bp-_E„) '-™,P,^'+1 "P" X e 

g-(£iv + i--E™)i V piEN + i+E^^2Ep)^__^ , 



p^m 



We find 



p^m 



The third contribution 



t:s = 



vl°\t,to), [Kit,to)/x'^\t,to) - Kito,to)\ vl^\t,to) 



'-Uta-ta) 



11 



is obtained similarly to 9tT2, permuting m and n. 
The fourth subleading contribution reads 

With some algebra we deduce 

tl (rors)>JV 



ii 

tl (r,s)>N 

tl 

+ E e"*''""''"^'*^*'' (Wn, V-W+l) (Mm, ^Ar+i)/^Ar+i,jv+ie-^" + ^ 
tl 

tl 



e 

tl 



fp f ^ ^ 



g-(B„-_E„)ti 



tl 



and we obtain 






„— (E'iV + l — -En) 1 

Cn,m.N+lhN+l,N+l 77^ ^r-. —6 " + ^ 

g— (-Ejv + i — -Em) ^ 

g— (-EjV + 1 — -Ernjt 



— C„ m,Ar+1^7V+l.W+l rp W~^ 7 

The last subleading contribution reads 

'/„(o) /^(2,i)/'+ ^„(n)\ /„(o) r^(2,i)Cf ^„(°) 



J 5 — t limn X 



( _r ., , ,p-(-Ejv+i--E™)to I „ -(£;jv+i--E„)to ] 



^Pm ^Pn 

t iljnn ^ 1 "T^ 

With to = t — 1 , we get 



We see that for n > m the dominating contribution Ti to eM^*'^ is in te ^«+i,"* with subleading terms 
T2 — Ts while ioT n < m the leading contribution is in g-(^N+i-s„)t^ 



12 

We have tested numerically our finding in the toy model of Ref. [19\, with ro-E„ = n, tq = 0.3, the 3x5 
matrix of couplings 

/0.92 0.03 -0.10 -0.01 -0.02^ 
i/, = (0|O,|n) = 0.84 0.40 0.03 -0.06 0.00 
\0.56 0.56 0.47 0.26 0.04 

and the hadronic matrix elements M„„ = 0.7— ^f, Mn n+m = „"" . The comparison between the analytical 
formulae and the numerical solution is plotted in Figure |8l It is encouraging to obtain such good agreement 
after i = 8. 
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FIG. 8: Analytical formulae for -R^rf ^^ compared to the numerical solution of our toy model. 
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